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■ Abstract 

O . 

The problem of constructing a necessary and sufficient condition for establishing 
the separability of continuous variable systems is revisited. Simon [R. Simon, Phys. 
Rev. Lett. 84, 2726 (2000)] pointed out that such a criterion may be constructed 
by drawing a parallel between the Peres' partial transpose criterion for finite di- 
mensional systems and partial time reversal transformation for continuous variable 
systems. We generalize the partial time reversal transformation to a partial scaling 
transformation and re-examine the problem using a tomographic description of the 
^ ■ continuous variable quantum system. The limits of applicability of the entanglement 

criteria obtained from partial scaling and partial time reversal are explored. 
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1 Introduction 



The qubit - an idealized two state quantum system - is the fundamental build- 
ing block out of which a universal quantum information processor may be en- 
visioned. Any representation of a qubit on a physical system is almost always 
an abstraction to the extent that we ignore all except two discrete states of the 
quantum entity. A complete description of the physical system representing a 
qubit usually involves continuous degrees of freedom like positions, momenta, 
relative phases and so on. It is fortunate that in most cases; especially in 
non-relativistic settings; we can ignore the continuous degrees of freedom of 
the qubit and focus on a subspace of its Hilbert space spanned by the eigen- 
states of an operator with a finite number of discrete eigenvalues like angular 
momentum, spin or polarization. 

Detecting, characterizing and understanding entanglement between finite di- 
mensional quantum systems like qubits and qutrits is an important piece in 
the developing paradigm of quantum information processing. Entanglement 
[1] is considered a physical resource that is used to construct quantum algo- 
rithms that can solve computational problems that classical algorithms and 
classical computers cannot realistically handle. 

Given a multipartite quantum system the question is whether it is possible 
to find a test that will detect entanglement that may be shared between the 
subsystems. Since separable states form a convex subset of the convex set of 
all possible states of the system, the Hahn-Banach theorem assures us that 
such a test must exist. Identifying the test is a non-trivial problem though. 
For two qubit systems Peres [2] suggested the partial transpose criterion which 
states that if a two qubit density matrix goes over to another density matrix 
(positive, trace one matrix) under partial transposition then the original state 
is separable. This was subsequently shown to be a necessary and sufficient 
condition for detecting entanglement in qubit-qutrit systems as well as two 
qubit systems by the Horodeckis [3]. A necessary and sufficient entanglement 
criterion for arbitrary multipartite quantum systems made up of finite dimen- 
sional subsystems is still not forthcoming in spite of all the interest that the 
problem attracts (see, for instance [4]). 

The Peres- Horodecki criterion is based on the fact that the transpose oper- 
ation is a positive but not completely positive map [5,6,7] on the state of a 
system. Choi [8,9,10,11] has shown that a not completely positive map applied 
to a part of an extended system will not be positivity preserving on all states 
of the whole system. For the combined state of two qubits, the transpose map 
applied to one of the qubits (partial transpose of the two qubit state) is not 
positive preserving on all entangled states. The same is true for the qubit- 
qutrit system. For two qutrit systems and higher, there exists certain states 
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called bound entangled states that remain positive under partial transpose 
even if its an entangled state and the Peres- Horodecki criterion ceases to be 
a sufficient condition for detecting entangled states. It is worth noting here 
that the transpose map applied to one of the subsystems is not a physically 
implementable operation. It is an algebraic manipulation that can be done 
on any given density matrix of the system to test for entanglement. So an 
alternate way of thinking about tests for entanglement is that we conceive of 
an operation (which may be unphysical) that makes sense only on separable 
states of the system and not on entangled ones. If an algebraic relation be- 
tween measurable quantities of the system that encapsulates the result of this 
operation can be found then that can be used to test entanglement. Several 
partial results pertaining to higher dimensional, multipartite state developed 
along similar lines can be found in the literature [12,13] 

The realization that the physical representation of an elementary qubit, even 
as something very simple like the spin of an isolated electron, must include at 
least one pair of canonical variables for a complete description is motivation 
to find entanglement criteria for continuous variable systems. Given that two 
qubits represented by, say, two electrons are in a "separable" state as far as 
the spin part of their combined wave function is concerned, one would like to 
ask the question whether they are really in a separable state given a complete 
description of the system in terms of spin, position and momentum coordi- 
nates. Of course, knowing the difficulties in finding such separability criteria 
even for generic finite dimensional, bipartite states, one might legitimately be 
suspicious of the chances of finding such a criterion for bipartite systems made 
of infinite dimensional subsystems. 

In [14], Simon pointed out that for bipartite continuous variable systems it 
is possible to find an analogue of the Peres-Horodecki criterion which is a 
necessary and sufficient test for entanglement in certain restricted cases. The 
partial transpose is shown to be equivalent, in the continuous variable case, 
to the time reversal of one of the subsystems. The partial time reversal is not 
only meaningful on separable states but in fact, in some cases, it is even imple- 
mentable using phase conjugation. On entangled states partial time reversal 
ceases to be a meaningful operation though. Simon finds an algebraic relation 
between the second moments of canonical position and momentum operators 
of the bipartite system that can be used to test for entanglement based on the 
partial time reversal transformation. 

In this paper we examine Simon's seminal result in greater detail and sug- 
gest a generalization based on partial scaling of the canonical variables. The 
partial time reversal is shown to be a special case of the partial scaling. We 
also phrase the discussion in terms of tomographic representation [15,16,17] 
of the state of the quantum system rather than using the Wigner function 
representation [18,19], the tomogram being the more readily measurable one 
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in an experimental setting. The scaling transformation can be applied to both 
the Wigner function and the symplectic tomogram [20,21,22] of the quantum 
state. The partial scaling of the tomogram has the advantage that it can be 
easily extended to states of systems with more than two pairs of continuous 
canonical variables. 

The test for entanglement devised using the partial time reversal and par- 
tial scaling transformations are based on the Heisenberg [23] and Robertson- 
Schrodinger [24,25] uncertainty relations. The connection of the Robertson- 
Schrodinger uncertainty relation to linear canonical transformations of the 
quadratures of quantum states with several canonical degrees of freedom is 
discussed in detail in [26,27,28]. Here we extend the discussion to the par- 
tial scaling transformations on arbitrary bipartite states and obtain a simpler 
derivation of Simon's algebraic relation that can be used as a test for entangle- 
ment. We also discuss the limits of applicability of the entanglement criteria 
based on both the partial time reversal and partial scaling. 

The organization of this paper is as follows: In section 2 we discuss the 
Robertson-Schrodinger uncertainty relations for multi-mode states. We in- 
troduce the partial scaling transformation on tomographic representation of 
quantum states in section 3. In section 4 we derive the criterion for entangle- 
ment based on the partial scaling transformation. The discussion in section 5 
is on the limits of applicability of the criterion and examples of systems on 
which the criterion will always work. Our conclusions are in section 6. 



2 Uncertainty relations for composite systems 

The quantum mechanical uncertainty relations expressed in terms of measur- 
able expectation values of combinations of the canonical observables associ- 
ated with a multi-mode system is the starting point for constructing a test 
for entanglement in such systems. We are interested in uncertainty relations 
that are invariant under canonical transformations of the system variables. 
The strategy is to then device a non-canonical transformation which, when 
applied to the relation, will indicate the presence of entanglement. 

Consider a system with N canonical degrees of freedom with N pairs of canon- 
ical variables denoted by £ a , a — 1, 2, . . . 2N. In our notation all £ a =2j (with 
even indices) are variables conjugate to £ a =2j-i- Canonical transformations 
are inhomogeneous, symplectic transformations on the variables of the form; 

£a ► S a p€/3 + c a- (1) 

These transformations are elements of the group Sp(2N, R) T(2N) where 
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Sp(2N, R) is the real symplectic group in 2N dimensions, T(2N) is the group 
of translations and represents the semi-direct product between the two. 
Canonical transformations leave the commutation relations between 6 invari- 
ant . 

For one degree of freedom, identifying the bilinear invariant under canonical 
transformation leads to the Robertson-Schrodinger uncertainty relation (units 
chosen so that h = 1): 

6 = 9 , 6 = p ; (6) = ci , (6) = c 2 

«& - <*)'>«6 - ^> 2 > - (mi±mi _ ClC2 y > i. (2) 

Note that the usual Heisenberg uncertainty relation 

((6- Cl ) 2 )((6-c 2 ) 2 )>^ (3) 



is not invariant under canonical transformations [28]. 

For simplicity in the discussion that follows we replace 6 — (£ Q ) with just 6 
and define Q = q-(q), P = p-(p), &qq = (Q 2 ), &pq = o QP = \(QP+PQ) = 
\({Q , P}) and dpp = (P 2 ). We can now write Eq. (2) in the form 



det C = det 



v + -a 



> 



(4) 



where 



V 



(6 2 ; 



1(66 + 66) 



V |(66 + 66) (6 2 ) 



yCTpQ <Jpp 



(5) 



and 



1 
-1 



(6) 



Rewriting the Robertson-Schrodinger uncertainty relation in terms of the dis- 
persion matrix V allows a direct generalization to systems with several canon- 
ical degrees of freedom. 



5 



For N degrees of freedom, 

V a f} = l({Z a ,Zf>}) ; a,(3 = l,2,...2N. 



(7) 



form a 2N x 2N matrix that transforms as an irreducible second rank tensor 
under the linear canonical (symplectic) transformations and have N invariants. 
The generalized uncertainty relations require that the matrix is non-negative 
and that each of its 2 x 2 minors be greater than or equal to 1/4. This condition 
can be written down as 

i 

C a p = V aj3 + 2^"/3 - (8) 



where £ is the canonically invariant block diagonal matrix, diag(f2, Q . . .Q). 
It is easy to see that 

det^>^. (9) 



While C a/ 3 and V a p are invariant under linear canonical transformations they 
are not invariant under scale changes on the £ a that are not contained in 
Sp(2iV, R). In particular under scaling C a p is not necessarily positive definite, 
much less satisfy the generalized Robertson-Schrodinger uncertainty bounds. 
This property can be used to construct a test for entanglement, starting from 
Eq. (8). But before looking at the test for entanglement we first look at how 
C a/ 3 can be computed for a given quantum state with continuous variables and 
how the scaling transformation may be implemented on it. Note that for one 
degree of freedom the scaled matrix V is always non-negative. The determinant 
of V being just a multiple of the determinant of the matrix before scaling. 
Therefore scaling will leave C positive as long as the overall multiplicative 
factor appearing before det V is greater than unity. 



3 Tomograms of quantum states and scaling transforms 



Given a density matrix p corresponding to a continuous variable system we 
can express it in the form of a distribution in phase space, p(p, q) where, 



Pip, q) 



q + 



X 



-ipx 



dx. 



Even thought p(p, q) appears to have continuous indices q and p, it is in fact 
only of discrete countably infinite dimension by virtue of the fact that p is in 
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a Hilbert space. So given p we can define a distribution in phase space. Such 
attempts like the Wigner-Moyal distribution given above and the diagonal 
coherent state representation (which is sometimes called the P-representation) 
have the problem that they are not guaranteed to be strictly non-negative. 
However there are two distribution functions that are strictly non-negative. 
The Husimi-Kano expectation values of p in a complete set of coherent states 
and the quantum tomogram [20,21]. The latter is defined by 

u(X,p,u) = tv[p5(X -pq- up)] (10) 



where the operator delta distribution may be understood in terms of the 
Fourier integral 

oo 

S(X -pq-up) = — ( ds jtix-w-vp) . (11) 
2ir J 



For systems with N degrees of freedom the definition of the tomogram may 
be generalized to 

N 

u(X, H, v) = tr[p Y[5(X S - p s q s - u s p s )\ (12) 

s=l 



The tomograms are all strictly non-negative being the integrals of the Wigner- 
Moyal distribution along the line pq + up = X. Note that even though ui 
appears to be the function of three variables, it is clear that the dependence 
on X can be written as 

u(X,p,u) = -!—uj(l,p/X, u/X). 



In our discussion we retain all three labels for reasons that will become clear 
later. 

The tomograms may be viewed as the line integrals of a density as in a clas- 
sical tomogram. Since they are always non-negative, the tomograms furnish a 
complete set of probability functions. It can be shown that the processes of 
computing the tomogram from a given p is invertible. If u(X,p, u) are given 
for all p and u we can find a linear formula for deducing the density distribu- 
tion and hence the density matrix. These are however rather complicated but 
still implementable (as in medical tomography). There is a direct inversion 
formula called the Bertrand-Bertrand formula [29,30]. 

For our present purposes we have to compute the second moments of canonical 
variables from the tomogram. By definition one has the tomographic dispersion 
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matrix elements, 

a XiXj {^u) = J dX(X t - (Xi)){Xj - (XJMX,!*,*) 



(13) 



where 



(Xi) = JdXXiUiiX, »,v). (14) 



From the tomographic dispersion matrix elements we obtain the elements of 
V a p using the following relations: 



a QiQi = a XjXj = 1 > = Vi = 0) 

G PjPj = a XjXj ( U 3 = 1 ' = Vi±j = 0) 

l r 

^QiPi = 2^^(^ = ^ = 1 ' AVj = *Vj = °) 

-(TXjXjilij = 1 , = Vi = 0) 

-crxjXjivj = l,fH = v m = 0)] 

OQ jPk = Xj X k {V>j = V k = \, fii^ij = V i+k = 0) (15) 



3.1 Scaling transform of the tomogram 



Now that we know how to compute the matrix elements of V from the tomo- 
gram we look at the effects of arbitrary scaling of the canonical variables on 
the tomogram. 

The scaling transform 

f(Qi,Pi) -»■ fs(qi,Pi) = J dq'dp' /C(?i,Pi;?-,p-)/(?-,p-) (16) 

with kernel 

N 

i=\ 

induces the following transformation on the tomogram: 

uj(X, fi,v)^uJs(X,n, u) 

r N 
= / dX' du' du' u(X', fx', v') JJ 6(Xi - xr) 

i=i 
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(18) 



The changes to the dispersion matrix elements brought about by the scaling 
transformation are given below: 



a QiPj- x x 

A Qi A Pj 

S ^QiQj 

a QiQj ~ \ \ 

s °PiP 3 / 1 Q\ 

a P i P j = ( 19 ) 
A Pi A Pj 



In equation (15) we computed elements of the dispersion matrix V by evalu- 
ating the elements of the tomographic dispersion matrix, crx-x-(f J, J u ) setting 
some of the /ij and equal to one and all others equal to zero. The elements of 
V after the scaling transformation is easily computed by setting the non-zero 
/ij and Vi equal to A^ 1 or A^. 1 (equal to the inverse scaling parameters of the 
corresponding canonical variables) rather than equal to unity. 

Given a tomogram uj(K, fi, u) we can compute the dispersion matrix V and use 
the inequality in Eq. (8) to test if the state corresponding to the tomogram 
satisfies the canonical uncertainty relations or not. If the tomogram is of a 
physical state we are assured that the inequality is satisfied. The same protocol 
can be employed to test whether scaled versions of the tomogram uj s (K, h, v) 
are acceptable to the extent that it produces dispersion matrices that respect 
the Robertson-Schrodinger uncertainty relations. 

In the case of a single mode, we know that scaling of the canonical variables 
does not change the positivity of V because the new matrix is just a positive 
multiple of the original one. In other words, det V — > A 2 detV. If u(X,fi,u) 
is such that the matrix C computed from it satisfies det C = det[V + > 
then, C computed from u s (X,fi,u) will also be such that detC > for all 
|A| > 1. It follows that if we have a separable multi-mode state given by the 
tomogram 

N 

u{X,n,u) =X> Jfc nw«(X i ,// i ,z/ i ) ; = ! (20) 

k i=l k 



then cu 5 (X, fi, v) obtained from arbitrary scalings of the canonical variables 
is always an acceptable tomogram provided cj(X, fi, v) leads to a dispersion 
matrix that satisfies (8). 
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Out of the 2N canonical variables {(q a ,Pa)}, & — 1, ---N associated with 
an N degree of freedom system, if we change the scale of some or all of the 
momentum variables to x a p a then that is not a canonical transformation. For 
separable states, the tomogram remains acceptable even under this partial 
scaling transform. For entangled states it turns out that this is not the case in 
general. In the next section we discuss how the partial scaling transform may 
be used to construct a necessary condition for the separability of multi-mode 
states. In restricted cases the condition is not only necessary but a sufficient 
test for entanglement as well. 



4 The criterion for separability based on partial scaling 

First let us look at the consequences of partial scaling on the dispersion matrix 
of a two mode quantum state. The state is described by the canonical variables 

= (qi, Pi, 22, P2), 

with 

Q% = q%- (qd ; Pi = Pi- (pd ; i = l, 2. 

We change the scale of the momentum variable of the second sub-system by 
a factor of x, i.e. 

P2 — ► — 21) 

x 

The dispersion matrix V2 of the second subsystem variables is transformed so 
that det V2 — > x 2 det Vi- So as long as \x\ > 1 we know that the partial scaling 
leaves the tomograms of separable two mode states acceptable (with positive 
C a fi). For a generic state the transformation changes the two mode dispersion 
matrix so that V s + |E is now 

/ „ „ _ \ 

2 



a QiQi + 2 a QiQ2 X( TQiP2 

°>iQi — 2 a P\P\ a PiQ-2 X(T PiP 2 

a Q 2 Qi a Q 2 Pi a Q 2 Q2 %< J Q 2 P 2 + I 

i 2 

\ x<j p 2 Qi x&p 2 p 1 xap 2 Q 2 — 2 x <Jp 2 p 2 ) 



(22) 



Let us denote the C and V matrices obtained after the partial scaling by 
C x and V x respectively, where x is the scaling parameter. To test for the 
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separability of a two mode state we test whether C x 
values of x with \x\ > 1. A consequence of V x + |E > is 



\/ x + f E > for all 



det 



V x + -E 



>0, 



which reduces to 

Ax 2 + 25a; + C>0 



(23) 



where 



A = det V - -(o-Q lQl ap lPl - Oq iPi ) 
B = ^( a QiP2 <J PiQ2 - a QiQ2 a PiP 2 ) 

C =Y^- ^Q2Q2 ( ^P2P2- (T Q 2 P 2 )- ( 24 ) 

For Eq. (23) to be always true, it is sufficient that the discriminant, 

B 2 - AAC < 0. (25) 



If we write the matrix in Eq. (22) as a block matrix, 



1 Vx V 12 ^ 



i 

+ 2 



n o 
o n 



(26) 



the condition (25) can be expressed as 

I det V 12 \ 2 - (4 det V - det Vi) (1-4 det V 2 ) < 0. 



(27) 



The inequality (27) must be satisfied by all separable two mode states and it 
need not be so for entangled ones. For certain classes of states (27) is satisfied 
only by separable ones and then the partial scaling transform furnishes a bona 
fide test for entanglement. We defer our discussion of this class of states to the 
next section and turn to the connection between partial scaling and Simon's 
[14] partial time reversal. 

The partial time reversal transformation for two mode states is a special case 
of the partial scaling with x = — 1. Time reversal of one of the subsystems of 
a bipartite quantum system is a canonical transformation as long as the state 
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of the system is separable. On entangled states the transformation ceases to 
be meaningfully defined. From (23), for separable states, if 



A + 2B + C > (28) 

then 

A - 2B + C > (29) 



also. Combining (28) and (29) we obtain 

11 1 

det V + — - -| det V 12 \ > - (det V x + det V 2 ) (30) 



Using 

det V = det V l det V 2 + | det V 12 \ 2 - tr[V l W l2 QV 2 QV&l] (31) 

we can re-write (30) in the form 

/l \ 2 
det Vi det V 2 + ( - - | det V 12 |J - tr^QV^QVaQV^fi] 

> i(det Vi + detV 2 ) (32) 

which is identical to Simon's criterion for separability of two mode states ob- 
tained using partial time reversal. The partial scaling transformation provides 
an alternate approach to deriving the criterion. 

Generalizing the criterion for separability based on the partial scaling trans- 
form to systems with N degrees of freedom is straightforward. Starting from 
the tomogram of a state of the system, ou (X, we first verify that it 

satisfies the Robertson-Schrodinger uncertainty relations by checking that 
Ca/3 = V a/ 3 + |S a/ 3 is a positive matrix. This may be done by making sure 
that all the principal minors of C are positive definite. We can now perform 
an arbitrary scaling described by the vector x = (x±, x 2 , . . .x 2 n) on the to- 
mogram. From the scaled tomogram we can now compute 

C^ = ^ + ^ 5 a,(3= 1,2,... 2N (33) 
where 

V* = [AcVAcW (34) 
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with _D X = diag(xi, X2, ■ ■ -X2n)- The 2N real parameters {x a } parameterize 
the Abelian scaling semi-group and we require that 

\xiX 2 \ > 1, \X 3 X A \ > 1, . . . , \x 2 N-IX2n\ > 1- 

The necessary condition for the separability of the state represented by the 
tomogram cu(X., fi, v) is that 

C x > (35) 
for all allowed choices of x. 

Out of the 2N scaling parameters we can always choose one pair, {xik-\i x ik) 
such that \xik~\Xk\ = 1 using the freedom to choose an overall scale factor 
that does not affect the positivity of C x . Furthermore, we are interested only 
in partial scalings that change the value of \x2j-iX2j\, j — 1, . . . , N. Without 
loss of generality we can scale only the momentum variables and leave the 
position variables unchanged so that all x 2 j-i = 1 and all x 2 j > 1. 

For two- mode systems, the choice x\ — X2 — £3 = 1, X4 = x^ 1 exhausts all 
the possibilities. For three mode systems there are more choices. By choosing 
X\ = £3 = £4 = £5 = Xq = 1 and x 2 = x~ l we can check whether the first 
mode is entangled to the remaining two. Each of the other two modes can 
similarly be tested for entanglement with the rest of the system. Tripartite 
entanglement can be tested for by scaling the momenta of two out of the 
three subsystems by choosing x\ = £3 = X5 = xq = 1 and X2 = x -1 , x 4 = 
The determinant of C x is now a polynomial in both x and y which is positive 
semi-definite in some domain in the xy-plane with \x\, \y\ > 1. 



5 Discussion 

The scaling transformation is not a canonical transformation and can be 
thought of as an effective scaling of the Planck's constant. A positivity of 
the density operator of a separable state is not sensitive to such scalings ap- 
plied to individual sub-systems. The entangled states, on the other hand, are 
sensitive to such scalings and in many cases this shows up by making C x 
negative for certain choices of x. 

Positivity of C x is a condition involving only the second moments of sym- 
metrized combinations of the canonical operators associated with the system. 
For distributions of continuous variables we require all the moments if we are 
to have complete information about the distribution. So positivity C x does 
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not guarantee that the distributions from which it came from corresponds to 
a separable state unless the second moments determine all the higher moments 
and the distribution itself. 



For symmetric distributions, all odd moments vanish and for a class of dis- 
tributions of which the most familiar is the Gaussian distribution, the second 
moment determines all the higher moments. So for systems characterized by 
Gaussian distributions the positivity of C x under partial scalings become a 
necessary and sufficient condition for separability. 

The construction of the matrix C x is related to the /—oscillator construction 
introduced in [31]. For /—oscillators the deformed annihilation operator has 
the form 



where and a are the standard creation and annihilation operators. If the 
function f(a^a) is a constant and equal to A -2 then the non-canonical trans- 
formation (36) corresponds to the Planck's constant being scaled to \h. Thus 
one can formulate the response of the density operator to scaling of the canon- 
ical variables in terms of the /—deformed oscillator with the function f(a)a) 
equal to a constant. 



5. 1 Limits of applicability of the criterion 

The criterion for separability can fail to detect entanglement in a mixture of 
two Gaussian distributions even if it is a necessary and sufficient condition on 
a single Gaussian. We construct an example of such a situation in this section, 
illustrating the limits of applicability of the partial scaling and partial time 
reversal transforms as the basis for constructing tests for entanglement. 

Consider a simple bipartite pure Gaussian state given by the wave function 



A = af(a t a) 



(36) 




1 



exp — -q T • M 1 ■ q 



(37) 




with 
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and 



\M\ = det M = m u m 22 - m 2 > 0. 

We have assumed for simplicity that all the elements of M are real. Since we 
assume for the purposes of the following analysis that we know the state, we 
may choose to describe it in terms of its tomogram, Wigner function or any 
other possible representation. The computation of the second moments of the 
distribution is transparent if we represent the state using its Wigner function, 
so we choose this option. The Wigner function corresponding to the Gaussian 
state is 



W(q,p) 



4tt 2 
1 

— exp 

TV 2 



ci 2 x #*(q + x/2)tf (q - x/2) e -* p ' : 



(q T • A!' 1 ■ q + p T • M ■ p) 



(3* 



Using the Wigner function we can compute the non-zero second moments of 
the state: 



(li) = ~7r~ , = ~Tr~ , (Q1Q2) = {mi) 



m 2 2 



m 
~2 



(39) 



and 



(P?> = 



2\M\ 



, (pi) 



2\M\ 



, (P1P2) = (P2P1) = 



m 



2\M\ 



(40) 



The Robertson-Schrodinger uncertainty relation for the state is given by 
/,„ .■ „ n \ 



£7=1 

2 





% 


m 





—i 


m.22 





m 




\M\ 


\M\ 


m 







i 





m 


—i 




\M\ 




|M| 



> 



(41) 



The pure Gaussian state in (37) is a minimum uncertainty state and we find 
that detC = as expected. Partial scaling, p 2 — > x~ l p 2 results in 



detC x = -(x-1) 



m 



16IMI 



< 



(42) 



which shows that the state is indeed entangled. 
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Now consider another pure Gaussian state 



* 2 (q) 



1 



cxp 



nJ\N 



--q T • N' 1 ■ q 



(43) 



with 



N 



n n 22 



, nn, n 2 2, n E R. 



A mixed state constructed out of the pure states in (37) and (43) of the form, 
p = ap 1 + (l- a)p 2 ; Pl = , < a < 1 (44) 



also satisfies the Robertson-Schrodinger inequality. The Wigner function for 
the mixed state reads, 



a 



W(q, p) = — exp 



7T^ 



m22 2 m n 2 2 2 

" W\ qi ~ W\ q2 ~ mnPl ~ m22P2 



m 



+2 IMI qiq2 ~ 2mpiP2 



1 — a 



IT' 



■ exp 



U 2 2 2 n H 2 2 2 

■7^1 ?1 - 777T 92 - n llPl ~ U 22P2 



n 



+2y77| 9i92 - 2npip 2 



(45) 



which, by inspection, is a non-separable function of the two pairs of canonical 
variables. The Robertson-Schrodinger uncertainty relations for the state is 
given by the inequality C m i X > where C m i X is the matrix: 



^ amn + (1 — a)nii i am + (1 — a)n 

• ttm22 _i_ (1— a)"22 n am (l-a)n 

1 \M\ ^ \N\ U 

am + (1 — a)n 

--^T TT— 



|M| |JV| 

am 22 + (1 - a)n 22 i 



am _ (l-a)n 
■|M| |iV| 



amn _i_ (1— a)nn 
|M| ^ |AT| / 



We find that 



det C m i x — 



a 2 (l-a) 2 (det[M - N]f 
16|M||iV| 



(46) 
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which is positive. Finding the determinant of the matrix C^ ix after an arbi- 
trary scaling p 2 — > x~ 1 p 2 has been applied is arduous and hardly enlightening 
but in the special case where x — — 1 when the partial scaling reduces to the 
partial time reversal we obtain 



detC* 



det C m i x 



am + (1 — a)n 



am (1 — a)n 



(47) 



We see that as long as 



am + (1 — a)n 



am (1 — a)n 



M 



N 



^ det Crnix 



the mixture of two entangled states appear to be separable with respect to 
the partial time reversal criterion. To view this as a failure of the test for 
entanglement one needs alternate ways of showing that the mixed state for 
which the Robertson-Schrodinger uncertainty relations are satisfied even after 
partial time reversal is in fact an entangled state. Since no such alternate 
method of showing that the state is entangled is forthcoming (otherwise we 
could just as well have used that as a test for entanglement) we look at a 
particular choice of M, N and a for which the test seems to fail and see what 
we may conclude from it. 

The choice mu = ran, m 22 = n 22 and m = —n (so that \M\ = \N\ but 
M ^ N) simplifies Eq. (47) to 

detC m „- — (48) 



Setting det C* nix = we obtain a fourth order equation in the parameter a. 
If we can find at least one solution a\ in the open interval (0, 1) we know 
that for a finite range of values of a between zero and one the scaled second 
moment matrix is such that the Robertson-Schrodinger uncertainty relations 
are satisfied. The equation det C^ ix = can be solved analytically using the 
change of variables a — > | — (3 reducing it to an equation in even powers of (5 
with solutions 



1 / 2|M| ± yl\M\{\M\+m?) 



4 \ m 



A plot of det C^ ix as a function of a for a particular choice of M and iV is 
given in Fig. 1. 
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a 



Fig. 1. det C!jni x for x = — 1 with m\\ = ri\\ = 77122 = "-22 = 0.5 and m = —n = 0.4. 

From Eq. (45) we know that for generic values of a the state is not separable. 
On the other hand from Fig. 1 we see that for a continuous range of values 
of a the partial time reversal (and partial scaling) criterion suggests that the 
state is separable. In the example we considered, by choosing m = —n, we are 
engineering the state so that the contribution to the cross correlations (q\q2) 
and (P1P2) from each of the two pure components of the state cancel each other. 
This is the reason why our test for entanglement based only of the second 
moments of the distribution fails. For generic values of a the mixed state we 
considered is entangled by construction. Since the sum of two Gaussian states 
is not another Gaussian we do not expect the second moments to determine 
all higher moments. We therefore expect signatures of the entanglement to 
show up in the nature of higher moments of the distribution. 



6 Conclusion 



We have derived a necessary (and in some cases sufficient) condition for the 
separability of quantum states of a multi-mode system with continuous vari- 
ables. The criterion is based on the partial scaling transform of the canonical 
variables of the system. We show how the scalings may be implemented on 
tomographic descriptions of the states. The test for separability is formu- 
lated in terms of the positivity conditions on the matrix of second moments 
computed from the tomogram. These conditions are based on the Robertson- 
Schrodinger uncertainty relations that all physical states must satisfy. Under 
partial scaling the matrix of second moments behaves differently for separa- 
ble and entangled states because partial scaling transforms a state to another 
physical state only if it is separable. The criterion for separability obtained 
using the partial scaling transformation reduces to the familiar result due to 
Simon [14] for separability of bipartite states when the scaling corresponds to 
a partial time reversal. 

For those states represented by distribution functions like the Gaussian for 
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which the second moment determines all the features of the distribution, the 
criterion derived here is a necessary and sufficient test for separability. We 
have also shown that the test is limited in its applicability because of the fact 
that it is based solely on the second moments and therefore cannot depend on 
all the features of the overall state. 
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